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Abstract. Using the effective potential approach for composite operators, we have 
formulated a general method of calculation of the truly non-perturbative Yang-Mills 
vacuum energy density (this is, by definition, the bag constant apart from the sign). It 
is the main dynamical characteristic of the QCD ground state. Our method allows one 
to make it free of the perturbative contributions ('contaminations'), by construction. 
We also perform an actual numerical calculation of the bag constant for the confining 
effective charge. Its choice uniquely defines the bag constant, which becomes free of 
all the types of the perturbative contributions now, as well as possessing many other 
desirable properties as colorless, gauge independence, etc. Using further the trace 
anomaly relation, we develop a general formalism which makes it possible to relate the 
bag constant to the gluon condensate defined at the same (3 function (or, equivalently, 
effective charge) which has been chosen for the calculation of the bag constant itself. 
Our numerical result for it shows a good agreement with other phenomenological 
estimates of the gluon condensate. We have argued that the calculated bag constant 
may contribute to the dark energy density. Its contribution is by 10 orders of magnitude 
better than the estimate from the Higgs field's contribution. We also propose to 
consider the bag energy as a possible amount of energy which can be released from the 
QCD ground state by a single cycle. The QCD ground state is shown to be an infinite 
and hence a permanent reservoir of energy. 
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1. Introduction 

In order to calculate physical observables from first principles in Quantum 
Chromodynamics (QCD) pQ it is not enough to know its Lagrangian. It is also necessary 
and important to know the true structure of its ground state. It is the response of 
the QCD vacuum which substantially modifies all the QCD Green's functions from 
their free counterparts. These full ("dressed") Green's functions are needed for the 
above-mentioned calculations. The vacuum of QCD is a very complicated confining 
medium and its dynamical and topological complexity means that its structure can be 
organized at various levels: classical and quantum [TJ |2j El SI El E] (and references 
therein). It is mainly non-perturbative (NP) by origin, character and magnitude, 
since the corresponding coupling constant is large. However, the virtual gluon field 
configurations and excitations of the perturbative (PT) origin, character and magnitude, 
due to asymptotic freedom (AF) pQ, are also present there. 

One of the main dynamical characteristics of the QCD ground state is the bag 
constant. Its name comes from the famous bag models for hadrons [7J[8], but its present 
understanding (and thus modern definition) not connecting to hadron properties. It is 
defined as the difference between the PT and the NP vacuum energy densities (VEDs) 
QUI QU [LJ. So, we can symbolically put B = VED PT - VED, where VED 
is the NP but 'contaminated' by the PT contributions (i.e., this is a full VED like 
the full gluon propagator, see below). At the same time, we can continue as follows: 
B = VED PT -VED = VED PT -[VED-VED PT + VED PT ] = V ED PT — [VED TNP + 
VED PT ] = —VED TNP > 0, since the VED is always negative. The bag constant is 
nothing but the truly NP (TNP) VED, apart from the sign, by definition, and thus is 
free of the PT contributions ('contaminations'). The symbolic subtraction presented 
here includes the subtraction at the fundamental gluon level, and two others at the 
hadronic level, i.e., when the gluon degrees of freedom should be integrated out (see 
section [3] below). In order to consider it also as a physical characteristic of the QCD 
ground state, the bag constant correctly calculated should satisfy some other necessary 
requirements such as colorlessness, finiteness, gauge-independence, no imaginary part 
(stable vacuum), etc. 

The main purpose of this paper is to formulate a formalism how to calculate 
correctly the quantum part of the bag constant, using the effective potential approach 
for composite operators [131 El H5]- in particular, to show how the above-mentioned 
subtractions are to be analytically made. On account of the confining effective 
charge, the bag constant has been numerically evaluated, satisfying all the necessary 
requirements mentioned above. Using further the trace anomaly relation [T6l fT71 fT8l fl9] . 
we also develop a general formalism which makes it possible to relate the bag constant 
to another important NP characteristic of the QCD ground state - the gluon condensate 
[TT] . Here, we do not use the weak coupling solution for the corresponding (3 function. 
Finally we present our numerical result for the bag constant, which is in a good 
agreement with other phenomenological estimates of the gluon condensate [Til [20] . 
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2. The VED 

The quantum part of the VED is determined by the effective potential approach for 
composite operators [131 [HI [15] • In the absence of external sources the effective potential 
is nothing but the VED. It is given in the form of the skeleton loop expansion containing 
all the types of the QCD full propagators and vertices, see Fig. [TJ So each vacuum 
skeleton loop itself is a sum of an infinite number of the corresponding PT vacuum 
loops (i.e., containing the point-like vertices and free propagators, see Fig. [2j where one 
term in each lower order is shown, for simplicity). The number of the vacuum skeleton 
loops goes with the power of the Planck constant, h. 




Figure 1. The skeleton loop expansion for the effective potential. The wavy 
lines describe the full gluon propagators D. The solid lines describe the full quark 
propagators S. T is the full quark-gluon vertex, while T3 and T4 are the full three- and 
four-gluon vertices, respectively. 




Figure 2. Infinite series for the gluon part of the VED (taking the first skeleton 
diagram in Fig. [1]). 

Here we are going to formulate a general method of numerical calculation of the 
quantum part of the TNP Yang-Mills (YM) VED in the covariant gauge QCD. The gluon 
part of the VED to leading order (the so-called log-loop level ~ h, the first skeleton loop 
diagram in Fig. [U and which PT expansion is shown explicitly in Fig. [2]) is analytically 
given by the effective potential for composite operators as follows [13] : 

V{D) = Uw) iTl {HD ° lD) " {D ° lD) + 1} ' (1) 

where D(q) is the full gluon propagator and Do(q) is its free counterpart (see below). 
The traces over space-time and color group indices are assumed. Evidently, the effective 
potential is normalized to V(D ) = 0, i.e., the free PT vacuum is normalized to zero, as 
usual. Next-to-leading and higher contributions (two and more vacuum skeleton loops) 
are suppressed at least by one order of magnitude in powers of h. They generate very 
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small numerical corrections to the log-loop terms, and thus are not important for the 
numerical calculation of the bag constant to leading order. 

The two-point Green's function, describing the full gluon propagator, is 

D^(q) = -i {T,M)d{-q 2 ^)+^L, p {q)} j 2 , (2) 

where d(—q 2 ;C,) is the gluon invariant function (dimensionless), the so-called Lorentz 
structure (sometimes, we will call it as the full gluon form factor or, equivalently, the 
effective charge ("running"), see below), while £ is the gauge-fixing parameter and 

Tf,u(q) = g^u - = 9^u ~ L^(q). (3) 

Its free PT counterpart D = D^ u (q) is obtained by putting the full gluon form factor 
d(— q 2 ;C) i n Eq. ([2]) simply to one, i.e., 

D%(q) = -i{T^(q)+£L^(q)}± (4) 

In order to evaluate the effective potential ([I]), on account of Eq. fl2]), we use the 
well-known expression 

Trlnp^D) = 8 x 41ndet( J D - 1 £>) = 32 ln[(3/4)d(-g 2 ; f) + (1/4)], (5) 

which becomes zero indeed when setting d(—q 2 ;^) = 1. 

Going over to four-dimensional Euclidean space in Eq. one obtains (e s = V(D)) 



e g = -16 



d 4 q 



ln[l + 3% 2 ;0]-^(g 2 ;0 + a 



(6) 



(2tt 

where constant a = (3/4) — 2 In 2 = —0.6363 and the integration from zero to infinity 
over q 2 is assumed. The VED e g derived in Eq. (jSJ) is already a colorless quantity, 
since it has been summed over color indices. Also it does not depend explicitly on the 
unphysical (longitudinal) part of the full gluon propagator due to the product (D^D), 
which, in its turn, comes from the above-mentioned normalization to zero. Thus it is 
worth emphasizing that the transversal ("physical") degrees of freedom only of gauge 
bosons contribute to this equation. Note, in the effective potential approach to leading 
order there is no need for ghost degrees of freedom from the very beginning in order to 
cancel the longitudinal ("unphysical") component of the full gluon propagator. This role 
is played by the normalization condition (that is why the ghost skeleton loops are not 
shown in Fig. [TJ. Furthermore, overall numerical factor 1/2 has been introduced into 
Eq. (1) in order to make the gluon degrees of freedom to be equal 32/2 = 16 = 8x2, 
where 8 color of gluons times 2 helicity (transversal) degrees of freedom, see Eqs. <^ 
and ©. 

In the connection with the above-mentioned normalization condition a few remarks 
are in order. It does not work for the higher order vacuum loops. As explained in Ref. 
|13j . for consistency with them in the PT QCD Green's functions, for example in the 
Hartree-Fock approximation, the Landau gauge should be used. In Ref. [2T] the effective 
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potential has been used to the two-loop order for the investigation of QCD chiral- 
symmetry breaking just in the Landau gauge and in the Hartree-Fock approximatiion. 
In the general case (i.e., beyond the PT and at any gauge), however, the cancelation 
of unphysical gluon modes should proceed with the help of ghosts as it is described in 
more detail in appendix A. 

The derived expression ([6]) remains rather formal, since it suffers from the two 
serious problems: the coefficient of the transversal Lorentz structure d(q 2 ; £) may still 
depend explicitly on £. Furthermore, it is divergent at least as the fourth power 
of the ultraviolet (UV) cutoff, and therefore suffers from different types of the PT 
contributions. 

3. The TNP VED 

In order to define the VED free of the above-mentioned PT contributions 
('contaminations'), let us make first the subtraction at the fundamental gluon level, 
namely 

d(q 2 ; = d(q 2 ; - d PT (q 2 ; + d PT (q 2 ; = d TNP (q 2 ) + d PT (q 2 ; 0, (7) 

where d PT (q 2 ; £) correctly describes the PT structure of the full effective charge d(q 2 ; £), 
including its behavior in the UV limit (AF, [TJ), otherwise remaining arbitrary. On the 
other hand, d TNP (q 2 ) defined by the above-made subtraction, is assumed to reproduce 
correctly the TNP structure of the full effective charge, including its asymptotic in 
the deep infrared (IR) limit. This underlines the strong intrinsic influence of the 
IR properties of the theory on its TNP dynamics. Evidently, both terms are valid 
in the whole energy /momentum range, i.e, they are not asymptotics. Let us also 
emphasize the principle difference between d(q 2 ;!;) and d TNP (q 2 ). The former is the 
NP quantity "contaminated" by the PT contributions, while the latter one being also 
NP, nevertheless, is free of them. Thus the formal separation between the TNP effective 
charge d TNP (q 2 ) and its PT counterpart d PT (q 2 ;£) is achieved. For example, if the full 
effective charge explicitly depends on the scale responsible for the TNP dynamics in 
QCD, say A 2 - the so-called mass gap (see section [5] below), then one can define the 
subtraction as follows: d TNP {q 2 ; A 2 ) = d{q 2 ; A 2 ) - d{q 2 ; A 2 = 0) = d{q 2 ; A 2 ) - d PT {q 2 ), 
which is, obviously, equivalent to the decomposition (jjj). In this way the separation 
between the TNP effective charge and its PT counterpart becomes exact, but not unique. 
Let us emphasize that the dependence of the full effective charge d(q 2 , A 2 ) on A 2 can 
be only regular. Otherwise it is impossible to assign to it the above-mentioned physical 
meaning, since A 2 can be only zero (the formal PT limit) or finite, i.e., it cannot be 
infinitely large. In principle, in some special models of the QCD vacuum, such as the 
Abelian Higgs model [221 E3] , the NP scale is to be identified with the mass of the dual 
gauge boson. 

There is also another serious reason for the subtraction in Eq. (JTj). The problem 
is that the above-mentioned UV asymptotic of the full effective charge may depend 
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on the gauge-fixing parameter £ explicitly, namely to leading order d(q 2 ,£) ~ (? 2^ 00 
( ln(g 2 /Ag C1) )) CO//bo , where the exponent (c /b ) < explicitly depends on the gauge- 
fixing parameter £ via the coefficient cq based on Ref. PQ, and Aq CD is the QCD 
asymptotic scale parameter. In this connection let us note that AF being a physical 
phenomenon does not depend on the gauge choice (it takes place at any gauge), while 
the UV asymptotic of the corresponding Green's function may be still gauge- dependent. 
This is just explicitly shown above. Evidently, in the decomposition (JTJ) just the PT 
part of the full effective charge will be responsible for this explicit dependence on the 
gauge choice. Subtracting it, we will be guaranteed that the remaining part will not 
depend explicitly on the gauge-fixing parameter (that is why the dependence on £ is 
not explicitly shown in d TNP (q 2 )). Let us note that if there is no exact criterion how 
to distinguish between the TNP and the PT parts in the full effective charge in Eq. (JTJ) 
as described above, then it is possible from the full effective charge to subtract its UV 
asymptotic only. However, in this case the separation between the TNP and the PT 
parts will be neither exact nor unique. For how to make this separation exact and 
unique at the same time see section |5j 

Substituting the decomposition (I7j) into Eq. (jHJ) and doing some simple 
rearrangements, one obtains 



7T 



dq 2 q 2 



ln[l + 3d TNP {q 2 )\ - ^ P (q 



(8) 



where the trivial integration over the angular variables in Eq. (6) has been already 
done. Here €pt is 

Sd PT (q 2 ; 



1 



epT 



71- 



dq 2 q 2 



In 1 



(9) 



1 + 3d TNP (q 2 

It contains the contribution which is mainly determined by the PT part of the full 
effective charge, d PT (q 2 ,£). The constant a should be also included, since it comes from 
the normalization of the free PT vacuum to zero. 

However, this is not the whole story yet. The first term in Eq. (8), depending 
only on the TNP effective charge, nevertheless, assumes the integration over the PT 
region up to infinity. It also represents the type of the PT contribution, which should 
be subtracted as well. If we separate the NP region from the PT one, by introducing 
the so-called effective scale q 2 ^ explicitly, then we get 



7T- 



q eff 



dq 2 q 2 



\n[l + 3d 1NP (q 2 )]-^d 



TNP I 2 



(l 2 



where 



~-pt 



71- 



-o ^ q 2 



eff 



ln[l + 3d TNP (q 2 )]~^d TNP (q 2 ) 



+e PT +e' PT , (10) 



(11) 



This integral represents the contribution to the VED which is determined by the TNP 
part of the full gluon propagator but integrated out over the PT region. Along with epx 
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given in Eq. it also represents a type of the PT contribution into the gluon part of 
the VED (jSJ), as mentioned above. This means that the two remaining terms in Eq. ffTOl) 
should be subtracted by introducing the TNP YM VED eym as follows: 

fFM = e g ~ e PT — CpTi (12) 

where the explicit expression for e Y M is given by the integral in Eq. (fID]) . 

Concluding, let us emphasize that both subtracted terms epx and e' PT , strictly 
speaking, are not the purely PT, since along with the nontrivial PT effective charge 
d PT (q 2 ) they contain the TNP effective charge d TNP (q 2 ) as well. So to call them the 
PT contributions is a convention. More precisely it is better to say that these terms are 
"contaminated" by the PT contributions. The above-mentioned necessary subtractions 
can be made in a more sophisticated way by introducing explicitly the ghost degrees of 
freedom (see appendix A). 

4. The bag constant 

The bag constant (the so-called bag pressure) is defined as the difference between the 
PT and the NP VEDs |3 HOI HH 112]. So in our notations for the YM fields, and as it 
follows from the definition by Eq. ( fl2i) . it is nothing but the TNP YM VED apart from 
the sign, i.e., 

Bym = —£ym = epr + e'pT ~ e g = 

i r q hf 

= - 2 dg 2 q 2 

7T 2 Jo 

where from now on we introduce the notation 

<f NP (q 2 ) - <*™ P (q*), (14) 

since d TNP (q 2 ) is the TNP effective charge aJ NP (q 2 ), as noted above. This is a general 
expression for any model effective charge in order to calculate the bag constant, or the 
TNP YM VED apart from the sign, from first principles. It is our definition of the 
TNP YM VED and thus of the bag constant. So it is defined as the special function 
of the TNP effective charge integrated out over the NP region (soft momentum region, 
< q 2 < Qeff)- ft i s f ree °f the PT contributions, by construction. In this connection, 
let us recall that e g is also NP, but 'contaminated' by the PT contributions, which just 
to be subtracted in order to get Eq. ffT3l) from Eq. (10). 

Comparing expressions (jHJ) and ( 1T31) . one comes to the following 'prescription 7 to 
get Eq. (flBl directly from Eq. (J6]): 

(i) Replacing d(q 2 ) — > d TNP (q 2 ) or equivalently, a s (q 2 ) — > a' FNP (q 2 ). 

(ii) Omitting the constant a which normalizes the free PT vacuum to zero. 

(iii) Introducing the effective scale q 2 ^ which separates the NP region from the PT one 
in the g 2 -momentum space. 

(iv) Omitting the minus sign for the bag constant. 



ln[l + 3a™ (q 2 )]-"-al» p (q 2 



(13) 
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At this stage the bag constant defined by Eq. (ITH1) is definitely colorless (color- 
singlet) and free of the PT contributions ("contaminations"). Let us remind that it 
also depends on only transversal degrees of freedom of gauge bosons (gluons). All 
its other properties mentioned above (finiteness, positivity, no imaginary part, etc.) 
depend on the chosen effective charge, more precisely on its TNP counterpart. It is 
worth emphasizing once more that in defining correctly the bag constant, three types 
of the corresponding subtractions have been introduced. The first one - in Eq. ([7j) at 
the fundamental gluon level and the two others - in Eq. ( TT2l) . when the gluon degrees of 
freedom were to be integrated out. 

For actual numerical calculations of the bag constant via the expression ( Tl3l) it is 
always convenient to factorize its scale dependence. For this purpose, let us introduce 
the dimensionless variable and the TNP effective charge as follows: 

a T s NP {q 2 ) = a T s NP {z) 1 where z = 4~ • ( 15 ) 

q eff 

From the general expression for the bag constant ffl3l) in these terms one then gets 

B YM{qlff) =<ltff X ^ym, (16) 

where we introduce the dimensionless TNP YM effective potential Qym, for convenience. 
Its explicit expression is 

•i 



&YM = — I d: : 



TV 2 







ln[l + SoF» p {z)\ ~ \<Z NP {z) 



;i7) 



Let us emphasize that in order to factorize the scale dependence in the effective potential 
it is necessary to choose the fixed scale, like of//; an< ^ n °t the scale which can be varied, 
for example like the mass gap which can go to zero in order to recover the PT limit (see 
section below). Eqs. f fl6|) and (|T7|) are the main subject of our consideration in what 
follows. It is worth emphasizing once more that these expressions are general ones in 
order to correctly calculate the Bag constant from first principles in any model gluon 
propagator. The only problem remaining to solve is to choose such TNP effective charge 
aJ NP (z) which, first of all should not explicitly depend on the gauge-fixing parameter 
£. At the same time, the implicit gauge dependence is not a problem. Such kind of 
the dependence is unavoidable in quantum or classical gauge theories, since the fields 
themselves are gauge-dependent PQ[2]- For the different TNP effective charges a TNP (z) 
one gets different numerical results. That is why the choice for its explicit expression 
(ansatz) should be physically and mathematically well justified (see below). 

In this connection, let us remind that the gluon Schwinger- Dyson (SD) equation 
is highly non-linear one, and it has a very complicated mathematical structure, so 
there is no hope for an exact solutions, the number of which is not even fixed 
PQ. This means that the number of independent solutions, obtained under specific 
truncation/approximation schemes and gauges, is not fixed a priori as well. From the 
very beginning they should be considered on equal footing. 
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5. Confining effective charge 

Let us choose the TNP effective charge as follows: 

af^V) -> af P (g 2 ) = f, (18) 

where the superscript "INP" stands for the intrinsically NP effective charge (for a such 
replacement see remarks below). Here A 2 = A 2 ^ is the so-called Jaffe-Witten (JW) 
mass gap, which is responsible for the large-scale structure of the QCD vacuum, and 
thus for its INP dynamics [21]. Let us note, that how the mass gap appears in QCD 
has been explicitly shown in our recent work in Ref. [25J. 

• The gauge independence is obvious, i.e., it does not depend explicitly on the gauge 
choice, since the mass gap is already renormalized, and hence it is a finite quantity. 

• It satisfies the Wilson criterion of confinement - area law for heavy quarks [26] [27] 
or, equivalently, leads to the linear rising potential between heavy quarks [281 ES] 
in continuous QCD, "seen" also by lattice QCD [301 |3_T]. In this connection 
a few remarks are in order. In the case of heavy quarks the response of the 
vacuum can be neglected, and therefore the interaction between them and gluons 
effectively becomes point-pike. Just this makes it possible to describe confinement 
of heavy quarks in terms of the linear rising potential, derived on the basis of 
the expression ( TT8|) . For the light quarks the response of the vacuum cannot be 
neglected. The corresponding quark-gluon vertex is not point-like, and therefore 
there is no way to analyze confinement of light quarks in terms of the linear rising 
potential. However, the expression ( 1TB]) can be still used for the solution of the SD 
equation for the quark propagator together with the corresponding Slavnov- Taylor 
(ST) identity for the vertex [32] • Confinement of light quarks is due to the analytical 
properties of the corresponding Green's functions (unlike the electron propagator, 
the quark propagator should have no imaginary part). This is a principle difference 
in the description of confinement for light and heavy quarks. 

• The functional dependence in the confining expression (1TB]) is, of course, the same 
for the YM fields and the full QCD. The dependence on the number of flavors can 
appear only in the mass gap. 

• It is exactly defined, since in the formal PT limit (A 2 = 0) the INP effective 
charge (1TB]) vanishes, and hence the bag constant itself. 

• It is uniquely defined as well. In order to show this explicitly, let us assume that it 
can be replaced by some arbitrary function as follows: 

«f(? 2 ;A 2 ) ^x/(g 2 ;A 2 ), (19) 

where f(q 2 ;A 2 ) is the dimensionless arbitrary function, which is regular at zero 
in order not to change confining properties of the INP effective charge (18). 
In this case it can be expand in Taylor series around small q 2 , so one obtains 
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f(q 2 ; A 2 ) = f{0) + {q 2 /M 2 )f'(q 2 ; A 2 ) + ..., where M 2 is some auxiliary mass squared. 
Then the INP effective charge in Eq. f|T9|) becomes 

ci N P (q 2 ; A 2 ) = ff(0) + ^f(q 2 ; A 2 ) + (20) 
and substituting this into the general decomposition ([7]), one finally obtains 

a s (q 2 ; A 2 ) = ai NP (q 2 ; A 2 ) + afV) = ^ + afV; A 2 ), (21) 

where not loosing generality we include the finite number f(0) into the mass gap, 
and retaining the same notation, for simplicity. The uniqueness is achieved at the 
expense of the PT effective charge, which now becomes regularly dependent on the 
mass gap (compare with the expression Qj). Evidently, the uniqueness is due to 
the singular at origin structure of the INP effective charge in Eq. ( IT8]) . In Ref. [33] 
it has been explicitly shown that the TNP part of the full gluon propagator as a 
function of the mass gap contains a regular at origin term as well. That is why 
it is not uniquely separated from the PT gluon propagator which effective charge 
is always regular at origin. We distinguish between the INP and the PT effective 
charges not only by the explicit presence of the mass gap, but by the character 
of the IR singularities as well [33] . So only after the replacement of Eq. (fl8l) the 
obtained expression for the bag constant ffTBl becomes free of all the types of the 
PT contributions ('contaminations'), indeed. 

• In our recent work [M] we have shown that the so-called INP gluon propagator is the 
purely transversal in a gauge invariant way, by construction. It exactly converges to 
the gluon propagator, which effective charge is in Eq. ffTB"]) . after the renormalization 
of the mass gap is completed. For preliminary analytical investigation of such 
behavior see Refs. [35l [36] as well (and references therein). Thus, we consider 
the expression (TT8l) not only as physically and mathematically well confirmed but 
as uniquely justified within the confining INP QCD [34J with its own mass gap 
identified with the JW mass gap for the pure YM fields (see above). 

• There also exist direct lattice evidences that the zero momentum modes are 
enhanced in the full gluon propagator (and hence in its effective charge) [371 E3 
39 . 40 1 (and references therein). A NP finite-size scaling technique was used in 
Ref. [JT] to study the evolution of the running coupling in the SU(3) YM lattice 
theory. At low energies it is shown to grow. The chosen analytical ansatz (18) can 
be considered as useful functional parametrization of these lattice results, indeed, 
while the scale of the enhancement is taken into account by the mass gap. 

• It is worth noting in advance that one of the attractive additional features of 
Eq. (18) is that it allows one to perform an analytical summation over the 
Matsubara frequencies in the generalization of the expression for the bag constant 
to non-zero temperatures. In this case one obtains the curve of the gluon matter 
pressure as a function of temperature. It and all other its derivatives (entropy 
and energy densities, etc.) then can be directly compared with the corresponding 
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thermal lattice QCD curves [121 H3]. This will make it possible for better 
understanding of the thermodynamical structure of the gluon matter (work in 
progress and preliminary numerical results are very encouraging). 

In conclusion, one may consider the expression ( JTB"]) as the confining ansatz, for 
simplicity. However, it is worth emphasizing that only it satisfies all the necessary 
conditions discussed above. Let us also note that for the theoretical and numerical 
results, depending on the confining effective charge, see discussion in section [HJ 



6. Analytical and numerical evaluation of the bag constant 

In terms of the variable in Eq. ( TT5l) for the INP effective charge ( ITBl) . one gets: 

z c q 2 A 2 

a I s NF {q 2 ) = a I s NP (z) = — , where z = — , and z c = — , 

Z q eff Qeff 

so that the dimensionless effective potential ffTTl) becomes, 

^ym(z c ) = — I dz z 
71 Jo 

Performing an almost trivial integration in this integral, one obtains 



Ml + (3z c /z)} - 



^ym{z c ) = Tn z l 



(23) 



(24) 



It is easy to see now that as a function of z c , the effective potential (124")) approaches 
zero from above as ~ z c at z c — > limit. At infinity z c — > oo it diverges as ~ — z c . 
At a fixed effective scale q 2 ^ and from Eq. ( 122]) it follows that z c — > is a correct PT 
regime, while z c — > oo is not a physical regime, since the mass gap A 2 is either finite or 
zero (the PT limit), i.e., it cannot be infinitely large. In other words, at a fixed effective 
scale one recovers the correct PT limit for the bag constant, i.e., the above-mentioned 
normalization condition is maintained for the bag constant, as it should be. 

The nontrivial second zero of the effective potential f l24j) follows obviously from the 
condition, 

3z c + 2 ln(l + 3z c ) - 18z 2 In (1 + (l/3z c )) = 0, (25) 

which numerical solution is 

z° c = 1.3786. (26) 

Evidently, through the relation f )22|) this value determines a possible upper bound for 
A 2 and lower bound for q^ff, since Bym/^ym is always positive/negative (see Figs. [3] 
and HJ). 

At z c = 0, i.e., A 2 = the effective potential ( l24j) vanishes identically, as it should 
be. From the above one can conclude that this effective potential as a function of z c has a 
maximum at some finite point, see Fig. |3j In the way how it has been introduced z c plays 
the role of the constant of integration of the effective potential though being formally 
a parameter of the theory. In general, by taking the first derivative of the effective 
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Figure 3. The ByM/Qeff effective potential vs. z c . The non-physical region is 
z c > Z®, since Bym should be always positive. At z c — the effective potential is also 
zero. 




Figure 4. The eym/^// effective potential vs. z c . The non-physical region is z c > z®, 
since €ym should be always negative. At z c = the effective potential is also zero. 



Vacuum Energy Density in the Quantum Yang - Mills Theory 



13 



potential with respect to the constant of integration one recovers the corresponding 
equations of motion [131 US]- Requiring thus 8VLym{.z c ) / dz c = 0, one obtains: 

^7 1 = 4 1n[l + (l/3^ c )], (27) 

which makes it possible to fix the constant of integration of the corresponding equation 
of motion at maximum. Its numerical solution is 

z rnax = q < ^ ( 2g ) 

so at maximum the ratio A 2 /q 2 *f is always less than one. At this point the numerical 
value of the effective potential ( 124|) is 

= 0.0263. (29) 

The bag constant defined in Eq. ffl6|) . and hence the corresponding INP VED ([TBI , as 
a function of q\^ or, equivalently, of the mass gap A 4 thus becomes, 

B Y m = ~e Y M = 0.0263 q$ ff = 0.1273 x A 4 , (30) 

where the relation 

q 2 eff = ( z ™ ax y l A 2 = 2.2 A 2 (31) 

has been already used. It is worth noting that a maximum for the bag constant 
corresponds to a minimum for the INP YM VED €ym (the so-called "stationary" state, 
see Fig. BJ. 

So, we have explicitly demonstrated that in the considered case the bag 
constant ( 130]) is finite, positive, and it has no imaginary part, indeed. It depends 
only on the mass gap responsible for the INP dynamics in the QCD ground state or, 
equivalently, on the effective scale squared separating the NP region from the PT one. 



O ( V 



max^ 



2vr 2 



max\ 



- In (1 + 3zl 



6.1. Scale-setting schemes and numerical results 

In order to complete the numerical calculation of the above defined bag constant all we 
need now is the value for the effective scale q 2 ft, which separates the NP region from 
the PT one. Similarly, the value for a scale at which the NP effects become important, 
that is the mass gap A 2 , also allows one to achieve the same goal. If the PT regime 
for gluons (as well as for quarks) starts conventionally from 1 GeV, then this number 
is a natural choice for the effective scale. It makes it also possible to directly compare 
our values with the values of many phenomenological parameters calculated just at this 
scale (see below). We consider this value as well justified and realistic upper limit for 
the effective scale defined above. Thus, using further the relation (j3T|) . one gets 

q 2 eff = 1 GeV 2 , and A 2 = 0.4564 GeV 2 . (32) 

Similarly, the numerical value of the mass gap A 2 has been obtained from the 
experimental value for the pion decay constant, F n = 93.3 MeV, by implementing a 
physically well-motivated scale-setting scheme [441 145] . In fact, we approximate the 
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pion decay constant in the chiral limit F® by its experimental value, since the difference 
between them can be a few MeV only. This is due to smallness of the corresponding 
light quark current masses. The pion decay constant is a good experimental number, 
since it is directly measured quantity, contrary to, for example the quark condensate or 
the dynamically generated quark mass. For the mass gap we have obtained the following 
numerical result A = 0.5784 GeV, so similarly to the relations fl32|) . one yields 

A 2 = 0.3345 GeV 2 , and q 2 eff = 0.733 GeV 2 . (33) 

In what follows we will consider this value as a realistic lower limit for the effective 
scale. One has to conclude that we have obtained rather close numerical results for the 
effective scale and the mass gap, by implementing rather different scale-setting schemes. 
It is worth emphasizing that the effective scale ( [33]) covers quite well not only the deep 
IR region but the substantial part of the intermediate one as well. 

For the above-mentioned possible upper bounds for A 2 and lower bounds for q 2 ^ 
our numerical results are for the scale-setting scheme (132]) : 

A 2 < 1.379 GeV 2 , and q 2 eff > 0.330 GeV 2 , (34) 

then similarly, based on the scale-setting scheme ( T33l : 

A 2 < 1.01 GeV 2 , and q 2 eff > 0.242 GeV 2 . (35) 

Evidently, their calculated values in each scale-setting scheme satisfy the corresponding 
bounds. 

For the bag constant (and hence for the INP YM VED) from Eq. (130]) . one obtains 

B YM = ~£ym = (0.0142 - 0.0263) GeV 4 , (36) 

where the first and second numbers in brackets correspond to the numerical values given 
in Eqs. ( 132]) and (!33|) . respectively. 

In conclusion, let us note that in the pure YM theory there is no way to calculate the 
mass gap independently of the well-motivated scale-setting scheme, that's the effective 
scale in this case, i.e., relations ( l32l) . The scale-setting scheme ( |33|) is based on the 
numerical value of the pion decay constant in the chiral limit. So this scheme is 
legitimated to use here as well, since the chiral quark condensates do not contribute to 
the VED in this limit, as it follows from the trace anomaly relation (see next section). 
For further discussion on the numerical value of Bym m different units see appendix B. 

7. The trace anomaly relation 

The TNP VED (and hence the bag constant) is important by itself as the main 
dynamical characteristic of the QCD ground state. Furthermore it assists in calculating 
such an important phenomenological parameter as the gluon condensate, introduced 
in the QCD sum rules approach to the physics of resonances The famous trace 
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anomaly relation [THl dZJ Q2J HH] in the general case of non-zero current quark masses 
is 

ew = ^r G % G % + E m ? w> ( 3? ) 

/ 

where 9 W is the trace of the energy-momentum tensor and G a being the gluon field 
strength tensor, while for the ratio (3(a s )/a s see discussion below. The trace anomaly 
relation which includes the anomalous dimension for the quark mass has been derived 
in Ref. [19], however, in our case of the pure gluon fields we can use the standard form 
of the trace anomaly relation (1371) . Sandwiching it between vacuum states and taking 
into account the obvious relation (0|9 MM |0) = 4e t , one obtains 

^ = (0|^G^G«J0) + 2>/<0|?/?/|0>. (38) 

/ 

Here et is the sum of all possible independent NP contributions to the VED (the total 
VED) and (0\qj-qf\0) is the chiral quark condensate. From this equation in the case of 
the pure YM fields (i.e., when the number of quark fields is zero Nf = 0), one can get 

(0|^G^|0)=4e yM , (39) 

where, evidently we saturate the total VED, e t by the TNP YM VED, eym defined in 
Eq. (TTBl) . i.e., putting e t = eyu + •••• Let us note that the same result, i.e., Eq. ( |39|) . 
will be obtained in the chiral limit for light quarks mj = 0, for / = 1, 2, 3 as well. 

If confinement happens then the [3 function is always in the domain of attraction 
(i.e., always negative) without IR stable fixed point pp. Therefore, it is convenient to 
introduce the general definition of the gluon condensate not using the weak coupling 
limit solution to the /3 function as follows: 

Thus, the above defined general gluon condensate will be always positive, as it should 
be. The importance of this relation is that it gives the value of the gluon condensate 
as a function of the bag constant whatever solution of the function in terms of a s 
is. However, let us remind that there is a correlation between the two sides of this 
equation. The bag constant, correctly defined in Eq. ( Fl3l) . depends, in general, on the 
TNP effective charge a™ p (q 2 ). On the other hand, the renormalization group equation 



(G 2 ) = -(0\^l G %G%\0) = 4 €ym — 4 Bym ■ (40) 



f ^f 1 = /?(«^ 2 )) (41) 



for the (3 function gives it in terms of the corresponding effective charge. This makes it 
possible to determine the ratio (f3(a s )/a s ) = (f3(a s (q 2 ))/a s (q 2 )), which appears in the 
left-hand-side of Eq. ( 140|) . Of course, this equation should be solved for the chosen TNP 
effective charge (see subsection 7.1). 
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Concluding, let us only note that the quantum part of the total TNP VED at 
log-loop level is: 

e t = e Y M + Nf e q , (42) 

where e q is the TNP quark skeleton loop contribution, see the corresponding skeleton 
loop diagram in Fig. [TJ It is an order of magnitude less than eyu because of much less 
quark degrees of freedom in the vacuum, and it is positive because of overall minus due 
to the quark loop. Evidently, in terms of the YM bag constant, one obtains 

e t = -B YM [1 - uN f ] , (43) 

where we introduce e q = vBym and i/<Cl. So the replacement of the total bag constant 
by its YM counterpart only is a rather good approximation from the numerical point 
of view. In this connection, let us remind that in the large iV c -limit the pure gluon 
contribution scales as N 2 , while the quark contribution scales only as N c [16]. However, 
in order to correctly calculate the bag constant in full QCD the quark part of the TNP 
VED e q is also important. Let us note that it is non-zero even in the chiral limit. This 
part will be investigated and calculated in our subsequent paper. 



(G 2 ) = -(0\^±G* u G° v \0) =AB YM = (0.0568 - 0.1052) GeV 4 . (44) 



7.1. Comparison with phenomenology 

Let us show explicitly now that our numerical values for the bag constant calculated 
in ( 1561 are in rather good agreement with the phenomenological values of the gluon 
condensate. Above we have already developed a general formalism which allows one 
to express the gluon condensate as a function of the bag constant. So substituting the 
numerical value of the bag constant into the Eq. (40), one obtains: 

i /3(a s 
4a s 

On the other hand, the renormalization group equation for the f3 function ( 14 ip after 
substitution of our solution for the INP effective charge (JIB"]) yields: 

(3(a s (q 2 )) = -a s (q 2 ), (45) 

as it is required for the confining theory where the /3 function should be always in 
the domain of attraction, i.e., negative (see in Ref. [T]). The corresponding ratio as it 
appears in the left-hand-side of Eq. ( 1441) is 

gW _ PMq 2 )) = _ L (46) 
a s a s (q 2 ) 

Substituting further this solution into the Eq. (144]) . it becomes 

(0|-G^G^|0) = AB YM = (0.0568 - 0.1052) GeV 4 , (47) 

which means that both sides of this relation between the Bag constant and the gluon 
condensate have been calculated by using the same expression for the INP effective 
charge, and hence for the corresponding function. So from the numerical point of 
view the Bag constant and the gluon condensate are in a self-consistent dependence from 
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each other, making thus the latter one free of all the types of the PT contributions. Our 
expression for the gluon condensate f|4"T|) allows one to recalculate any gluon condensate 
at any scale and any ratio, (3(a s )/a s . To the gluon condensate a physical meaning 
can be indeed assigned as the global (average) vacuum characteristic which measures a 
density of the TNP gluon fields configurations in the QCD vacuum. 

However, it cannot be directly compared with the phenomenological values for the 
standard gluon condensate estimated within different approaches [20]. The problem is 
that it is necessary to remember that any value at the scale as in Eq. ( 1331) (lower bound 
in the right-hand-side of Eq. ( I4"T|) ) is to be recalculated at the 1 GeV scale. Moreover, 
both values explicitly shown in Eq. (147p should be recalculated at the same ratio, as 
mentioned above. 

In phenomenology the standard ratio of the gluon condensate and its numerical 
value is: 

G 2 = (-G 2 ) = (0\^LG%G%\0) « 0.012 GeV 4 , (48) 

which can be changed within a factor of ~ 2 [TT] (let us recall that this ratio comes 
from the weak coupling solution for the (3 function, see for example in Ref. |47J). 

Thus in order to achieve the same ratio the both sides of Eq. (1471) should be 
multiplied by 4(a s /ir). For the numerical value of the strong fine structure constant 
we use a s = a s (rriz) = 0.1187 from the Particle Data Group |48j. In addition, the lower 
bound should be multiplied by the factor (1/0. 733) 2 = 1.86, coming form the numerical 
value by Eq. (1331 . Then the recalculated gluon condensate in (f47|) . which is denoted as 
G2, finally becomes (i.e., both numbers in Eq. (T4T1) coincides, as it should be) 

G 2 « 0.016 GeV 4 . (49) 

This numerical value for the gluon condensate should be compared with the numerical 
value coming from the phenomenology, see Eq. ( T481) above. This shows that all our 
numerical results are in good agreement with various phenomenological estimates [TTJ 
[20] . taking into account that the quark contributions are approximately an order of 
magnitude less than the pure YM one to the full bag constant (see remarks in this 
section just before subsection 7.1). This confirms that our numerical values for the bag 
constant and hence for the gluon condensate are rather realistic ones. 

8. Conclusions 

In summary, we have formulated a general method how to calculate numerically the 
quantum part of the TNP YM VED (the YM bag constant, apart from the sign, by 
definition) in the covariant gauge QCD from first principles. For this purpose we have 
used the effective potential approach for composite operators to leading order [13] . It has 
an advantage to be directly the VED (the pressure) in the absence of external sources. 
The bag constant is defined as a special function of the TNP effective charge integrated 
out over the NP region (soft momentum region), see Eq. ({TBI . At this stage the bag 
constant is colorless (color- singlet) and depends only on the transversal ("physical") 
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degrees of freedom of gauge bosons. It is also free of the PT contributions by its 
construction. This has been achieved due to the subtractions at the fundamental level 
as given by Eg. ([7]), as well as due to all other subtractions explicitly shown in Eq. ( TT2l) . 
when the gluon degrees of freedom were to be integrated out. Thus, our equations (TL6I) 
and ( TT71) are general ones in order to correctly calculate the bag constant as a function 
of any properly defined TNP effective charge. 

For the concrete calculation of the bag constant we replace the TNP effective 
charge by its confining INP counterpart in Eq. ( fl8l) . since it is exactly and uniquely 
separated from the PT effective charge. The INP effective charge depends regularly 
on the mass gap, which is responsible for the large-scale structure of the QCD ground 
state [2U ES]. The scale-setting schemes have been chosen by the two different ways, 
leading, nevertheless, to a rather close numerical results for the mass gap and hence 
for the effective scale. The calculated bag constant in addition, is: finite, positive, 
and it has no imaginary part (stable vacuum). It is also a manifestly gauge- invariant 
quantity (i.e., does not explicitly depend on the gauge-fixing parameter as it is required). 
The separation of "soft versus hard" gluon momenta is also exact because of the 
maximization/minimization procedure. It becomes possible since the effective potential 
( |24|) as a function of the constant of the integration z c and hence of the mass gap A 2 has 
a local maximum, see Fig. [3j This also makes it possible that in the above-mentioned 
scale-setting schemes either the mass gap or the effective scale is only independent, since 
the other one is to be determined via the relation f l3Tj) . In the scale-setting scheme (|32|) 
the effective scale is independent, while in the second scale-setting scheme ( |33i) the mass 
gap is independent. It is worth emphasizing that the bag constant in our approach is 
not simply the difference between the PT and NP VEDs, which is finite, colorless and 
manifestly gauge-invariant, etc. It is the energy density (apart from the sign) of the 
system of stable configurations of the purely transversal quantum virtual fields with the 
enhanced low-frequency components /large scale amplitudes due to the NL interaction of 
massless gluon modes, and which is being at "stationary state", i.e., being in the state 
with the minimum of energy, see Fig. [^J 

In order to compare our numerical results with phenomenology we develop a general 
formalism which makes it possible to relate the bag constant to the gluon condensate in 
a unique and self-consistent way. In other words, the gluon condensate is defined and 
calculated at the same effective charge, which has been chosen for the calculation of the 
bag constant. For this purpose we use the trace anomaly relation without applying to 
the weak coupling solution for the corresponding /3 function. In its turn, it is a solution 
of the corresponding renormalization group equation for the effective charge Eq. (41). 
Our numerical results turned out to be in good agreement with phenomenological 
values of the gluon condensate calculated and estimated within different approaches 
and methods [TTII20]. 

It is instructive to briefly summarize our theoretical and numerical results for the 
bag constant in general and our specific ways: 
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General properties of the bag constant determined by Eqs. (16)-(17) are: 

• colorless (color-singlet); 

• electrically neutral; 

• transversal, i.e., depending only on "physical" degrees of freedom of gauge bosons; 

• free of the PT contributions ('contaminations'). 

Results, depending on the confining effective charge Eq. (18) are: 

• the explicit gauge invariance; 

• uniqueness, i.e., it is free of all the types of the PT contributions now; 

• finiteness; 

• positiveness; 

• no imaginary part (stable vacuum); 

• existence of the stationary state for the corresponding YM energy density (negative 
pressure, see Fig. 4); 

• the final dependence on the mass gap only; 

• a good numerical agreement with phenomenology. 

The above remarkable features all together are unique. Apparently, it is due to 
the confining expression ( jl~8l) and the correct determination of the bag constant itself 
in this investigation. It has been made in accordance with its modern definition as the 
difference between the PT and the NP VEDs [9j M, EH 02] • 

Our method can be generalized on the multi-loop skeleton contributions to the 
effective potential approach for composite operators, as well as to take into account 
the quark degrees of freedom, as plotted in Fig. [TJ These terms, however, will produce 
numerical contributions an order of magnitude less, at least, in comparison with the 
leading log- loop level gluon term given by Eq. (CQ). What is necessary indeed, is to 
be able to extract the finite part of the TNP VED in a self-consistent and manifestly 
gauge-invariant ways. This is provided by our method which thus can be applied to any 
QCD vacuum quantum and classical models at any gauge (covariant or non-covariant). 
It may serve as a test of them, providing an exact criterion for the separation "stable 
versus unstable" vacua. Using our method we have already shown that the vacuum 
of classical dual Abelian Higgs model with string and without string contributions is 
unstable against quantum corrections [19], [50] . 

It would be also interesting to apply our general equations (fT6|) and ( |T7I) in order to 
calculate the bag constant within the recently obtained analytical results for the gluon 
propagator in Refs. [SH E2J [53] EH EE] [56], [57] [58] . The general formalism developed in 
our paper is aimed first of all at the analytical calculations of the bag constant (or the 
vacuum energy density) in any model gluon propagator in continuous QCD. However, 
as mentioned above the chosen ansatz (18) can be considered as useful parametrization 
of the corresponding lattice results. In this way our formalism can be extended to the 
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lattice calculations as well. Choosing an appropriate parametrization of any lattice 
result for the gluon propagator (there is a lot of recent lattice data [591 EOl EE E2] and 
references therein), one then can substitute it into our analytical expressions (16)-(17). 
Such a combination of the lattice and analytical calculations can be rather effective 
indeed, in order to understand what is the physics behind the lattice numbers and 
curves. On the other hand, all the analytical expressions and calculations will be put on 
solid numerical ground provided by the lattice simulations. So there is no doubt that 
the analytical and lattice calculations should not exclude each other, but contrary they 
should complement each other. All these possible developments are, of course, beyond 
the scope of the present investigation, and they have to be done elsewhere. 
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Appendix A. The general role of ghosts 

Let us begin with recalling that due to the above-mentioned normalization condition in 
the initial Eq. (1), its elaborated counterpart in Eq. (6) depends only on the transversal 
("physical") component of the full gluon propagator. So there is no need for ghosts to 
cancel its longitudinal (unphysical) component, indeed. However, it is instructive to 
discuss the role of ghosts in general, and to clearly show that their explicit introduction 
leads to the same result for the bag constant, in particular. 

Following Ref. [13], the effective potential at the same log- loop order for the ghost 
degrees of freedom analytically is: 

V(G) = -iJ -^TriHGo'G) - (G^G) + 1}, (A.l) 

where G = G(k) = i/k 2 (l + b(k 2 )) is the full ghost propagator, where b(k 2 ) is the 
ghost self-energy, while Go = Go(k) = i/k 2 is its free PT counterpart. Trace over color 
group indices is assumed. Evidently, the effective potential is normalized to V(G ) = 
in the same way as the gluon part in Eq. ([[]). Substituting these expressions into the 
ghost term (1A.1I) and again doing some algebra in four- dimensional Euclidean space, one 



formally obtains that V(G) = e g h = f dk 2 f{b{k 2 )). This, in general, divergent constant 
contribution should be of course, regularized in order to assign to it a mathematical 
meaning. So the explicit functional dependence of the ghost propagator /self-energy on 
its argument is not important, since within the effective potential approach to calculate 
the VED it is always only constant. We have to sum up all the contributions for the 
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pure YM fields at the same skeleton log-loop order. The relation given by Eq. ( II 01) then 
should look like as: 



1 r lff J 2 2 

e 9 + e 9 h = 2 / d(? q 

7T In 



Ml + 3d™ V)] - \<F NP {?) 



+ 

+ e PT + e' PT + e gh . (A.2) 

It is worth emphasizing that, the right-hand-side of this relation may still suffer from 
unphysical singularities by the integral in Eq (Q, defining epx- The problem is that 
the PT effective charge, d PT (q 2 ), which is responsible for AF in QCD at large q 2 (see, 
for example our paper [31]), may have, in general, unphysical singularities below the 
scale A.q CD , since in the integral ([9]) the integration is from zero to infinity. In addition, 
as mentioned above the integral pip , defining e' PT , may be still divergent. Thus the 
left-hand-side of the relation ( 1A.2I) is formal one, indeed. It suffers from various types 
of unphysical singularities which may appear in its right-hand-side. In order to get a 
physically meaningful expression, one has to remove the two integrals (Q and (TTTT) from 
Eq. (EI). This is to be done with the help of a ghost term by imposing the following 
condition of cancelation of unwanted terms in the most general form: epT + e' PT +e gh = 0. 
This condition can be always fulfilled, since it is a relation between three different 
(unknown in general) regularized constants. Then the relation (1A.2|) thus becomes: 



£yM — e g ~ e PT — t'pT 

in complete agreement with the relation ( fl2]) . and hence with the definition of the 
bag constant ffT3"j) . as it should be. So the TNP gluon contribution to the VED has 
been determined by subtracting unwanted terms by means of the ghost contribution. 
Evidently, the subtracted terms are of no importance, while a ghost term plays no 
explicit role for further consideration. 

In QCD the general role of ghost degrees of freedom is to cancel all the unphysical 
degrees of freedom of gauge bosons [TJ [63] , maintaining thus unitarity of the S- matrix. 
This is the main reason why they should be taken into account together with gluons 
always. This means that nothing should explicitly depend on them after the above- 
mentioned cancelation is performed. One of the main purposes of their introduction is 
to exclude the longitudinal (unphysical) component of the gluon propagator in every 
order of the PT, thus going beyond it and thus being a general one, indeed. If there 
is no need to cancel the longitudinal component of gauge boson propagators, then they 
should be used to eliminate the unphysical singularities of gauge bosons below the QCD 
asymptotic scale (as it was described above), or some other ones which may be inevitably 
present in any solution/ansatz for the full gluon propagator. If one knows the ghost 
propagator exactly, then the above-mentioned cancelation of unphysical singularities of 
gauge bosons should proceed automatically, as usual in the PT calculus (if, of course, 



ln[l + 3^^)]-^— (q' 



(A.3) 
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all calculations are correct). For such an exact cancelation of the longitudinal part 
of the gluon propagator by the free PT ghost propagator in lower order of the PT 
see, for example Ref. [63J. But if it is not known exactly or known approximately 
(depending on the truncation/approximation scheme), as usual in the NP calculus then 
nevertheless, one has to impose the corresponding condition of cancelation in order to 
fulfill their general role. This just has been done above. Thus our subtraction scheme 
is in agreement with the general interpretation of ghosts to cancel all the unphysical 
degrees of freedom of gauge bosons [U [63] . So by themselves the ghosts cannot change 
the truly NP dynamics of QCD, associated with the transversal component of the full 
gluon propagator in Eq. (2) and described by its Lorentz structure or, equivalently, by 
its effective charge (see Ref. [25] as well). 

Whatever solution(s) for the full gluon propagator obtained by lattice QCD 
[59| EUJ EU [62] (and references therein) and by the analytical approach based on the 
corresponding SD system of equations [5H [52J, [53J, EU [551 123 EO EH] (and references 
therein) might be (smooth, singular, massive, etc.), it, however, should not undermine 
the above-mentioned general job of ghosts. It is worth emphasizing that by no 
coincidence in all the papers cited above the transversal Landau gauge has been chosen 
by hand from the very beginning. So there is no and cannot be the explicit dependence on 
the ghost degrees of freedom in any expressions for the physical quantities, in general, 
and in the expression for the bag constant, in particular. In this connection, let us 
remind that the confining effective charge ([TBI is the effective charge of the relevant 
gluon propagator, which becomes the purely transversal in a gauge invariant way, by 
construction [33| 134]. 

Nevertheless, this does not mean that we need no ghosts at all. First of all, we need 
them in the higher orders of the two-particle irreducible vacuum graphs in the skeleton 
loop expansion of the effective potential [13], since for them the simple normalization 
of the free PT vacuum to zero does not work. So the cancelation of unphysical gluon 
modes should proceed with the help of the ghost degrees of freedom, as it was described 
in this appendix above. It is necessary to understand that the transversality of the 
gluon propagator in the Landau gauge in order to correctly treat the PT QCD Green's 
functions without ghosts is not enough to insure unitarity of the S'-matrix in QCD. The 
whole machinery of all the ST identities and the corresponding SD equations is still 
needed in order to insure the unitarity cancelations even in the Landau gauge. 

For example, the quark ST identity, contains the so-called ghost-quark scattering 
kernel explicitly [lj. This kernel still makes an important contribution to the identity 
even if the gluon propagator is transversal [641 [65] • Omitting ghosts at all in this 
identity, one will lose an important piece of information on the quark degrees of freedom 
themselves. As a result, any solution of the quark SD equation will suffer from unphysical 
singularities in the complex momentum plane. The problem is that via the quark-gluon 
vertex this equation will crucially depend on the term which comes from the identity even 
if the gluon propagator is transversal. The completely NP analysis of this identity on 
the basis of the double pole structure of the full gluon propagator in the IR, Eq. (18), has 
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been made in our earlier papers [SU EZ] • We have derived the corresponding expression 
for the quark-gluon vertex following Ref. [68J only in more sophisticated way (see Ref. 
[32] as well). We will take this result into account when we will directly calculate the 
confining quark contribution to the bag constant, as mentioned in section 7 just before 
subsection 7.1. 



Appendix B. Numerical values for Bym m different units 

In order to show explicitly what magnitude of numbers we are dealing with, let us 
present our numerical value for the bag constant given by Eq. ( 136]) in different units, 
namely: 

Bym = -gym = (0.0142 - 0.0263) GeV 4 
= (1.84 - 3.4) GeV/fm 3 

= (1.84 - 3.4) x 10 39 GeV/cm 3 . (B.l) 

This is a huge amount of energy stored in one cm 3 of the QCD vacuum even in "God- 
given" units H = c = 1. Using the number of different conversion factors (see, for 
example Ref. [63] or the particle data group [48]) the bag constant can be expressed in 
different systems of units (SI, CGS, etc.). 
Taking further into account that 

1 GeV = 1.6 x 10 10 J = 4.45 x 1(T 23 GWh, (B.2) 
from Eq. f lBlj) one finally gets (1 W = 10" 3 kW = 10" 6 MW = 10~ 9 GW) 

B YM = (8.2 - 15) x 10 16 GWh/cm 3 (B.3) 
or, equivalently, 

E YM = B YM cm 3 = (8.2 - 15) x 10 16 GWh ~ 10 17 GWh (B.4) 

in familiar units of watt-hour (Wh). Let us note that if one puts the effective scale 
squared as small as realistically possible q 2 e ^ = 0.242 GeV 2 (see Eq. (1351) ). then 
the previous number will be only slightly changed, namely Eym = Bym cm 3 = 
(4.8 — 8.7) x 10 15 GWh. So both numbers still indicate a huge amount of the bag 
energy E Y m stored in one cm 3 of the QCD vacuum. 

It is especially interesting to compare these numbers with the total production of 
primary energy of the 25 EU countries in year 2004 which was [69] (see also Ref. |70j ) 

E t ~ 10.2 PWh = 10.2 x 10 6 GWh ~ 10 7 GWh, (B.5) 

where 1 PWh = 1 Petawatt-hour. Approximately 1/3 of this energy was produced by 
nuclear power plants [691 EQ]- The huge difference between the numbers in Eqs. (IB. 41) 
and ( IB. 51) is very impressive and leads to some interesting still speculative but already 
possible discussion in appendix D below and in our preliminary work [71] . 
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Appendix C. Contribution of Bym to the dark energy problem 

Apparently, our bag constant ( IB. II) may also contribute to the so-called dark energy 
density [72J. At least, from the qualitative point of view it satisfies almost all the criteria 
necessary for the dark energy/matter candidate (see here section 8 and discussions in 
Refs. [721 [73]). From the quantitative numerical point of view it is also much better than 
the estimate from the Higgs field's contribution to the VED, which is about [741 175] 

q h ~ 10 8 GeV 4 . (C.l) 
In this notation our value (B.l) is about 

Oour ~ 1(T 2 GeV 4 . (C.2) 
The observed VED is very small indeed, namely 

£W ~ lO" 46 GeV 4 , (C.3) 

see Refs. [7H [751 [76] . So relatively to the value inferred from the cosmological constant 
(i.e., the above-mentioned observed VED) 

Qui Qvac ~ 10 54 , (C.4) 

while our is 

Qour I Qvac ~ 10 44 , (C.5) 

i.e, some 10 orders of magnitude better, which is expected from the direct comparison 
of the estimate (C.l) with our value (C.2). 

Let us note that calculating at the Plank length scale [48] , we will obtain the same 
ratio, as it should be. From Eq. (B.l) it follows that 

g our ~ 10 39 GeV/cm 3 = 1(T 60 GeV/L 3 , (C.6) 

where we used cm = 10 33 L p and L p denotes the above-mentioned Plank length [48J. In 
this units the observed VED is 

g mc ~ 10" 46 GeV 4 ~ 1(T 5 GeV/cm 3 = 10~ 104 GeV/L 3 , (C7) 

so that the ratio between (C.6) and (C.7) becomes again (C.5), indeed. Of course, 
the ratio (C.5) still remains very large, but it is much better than the ratio (C.4), as 
emphasized above. Other possibility how QCD can be related to the dark energy puzzle 
has been described in Ref. [77] (and references therein). 

Concluding, the vacuum for which the value (C.3) has been measured should not 
be mixed with the vacuum of any quantum field gauge theory. For the former one its 
energy is always positive (i.e., above zero), so the vacuum is simply treated as empty 
space. The energy of the latter one is always negative (i.e., below zero), and it is full 
of any kind of quantum excitations, fluctuations, etc. However, the QCD bag constant 
is always positive, finite, gauge- invariant, etc. (if it has been correctly defined and 
calculated like in this work) . That is the primary reason why we can compare our value 
(C.2) and the estimate (C.l) with (C.3). 
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Appendix D. Energy from the QCD vacuum 

The Lamb shift and the Casimir effect are probably the two most famous experimental 
evidences of zero-point energy fluctuations in the vacuum of Quantum Electrodynamics 
(QED) [TU EE EQj IH]. Both effects are rather weak, since the QED vacuum is mainly 
PT by origin, character and magnitude (the corresponding fine structure constant 
is weak). However, even in this case attempts have been already made to exploit 
the Casimir effect in order to "observe" the negative energy and related affects [80] 
and even to release energy from the vacuum (see, for example Refs. [821 [83] and 
references in the above-mentioned reviews [781 179]). In Ref. [HI] by investigating 
the thermodynamical properties of the quantum vacuum it has been concluded that 
no energy can be extracted cyclically from the vacuum (see, however Ref. [78] and 
references therein). Let us also note that in QED the photon propagator always remains 
PT even "dressed" [251 Ell ESI EH]. So formally we can define the bag constant in this 
theory as B QED = VED - VED = -VED > 0, since VED = VED(D ) = in the 
effective potential approach to leading order [13]. It would be interesting to perform 
such a calculation, which will give one a correct finite value of the VED in QED, if, or 
course, the above proposed definition of the QED bag constant makes sense. But it is 
beyond the scope of the present investigation, and should be done elsewhere. 

Since the QCD fine structure constant is strong, the idea to exploit the QCD 
vacuum in order to extract energy from it seems to be more attractive. However, before 
discussing the ways how to extract, it is necessary to discuss which minimum/maximum 
amount of energy at all can be released in a single cycle. Who thinks that it is too early 
to discuss such kind of topic (though we do not think so) may entirely skip this appendix. 

The bag constant calculated here is a manifestly gauge- invariant, real and colorless 
(color- singlet) quantity, i.e., it can be considered as a physical quantity. In fact, in 
this paper we have formulated a renormalization program to make the bag constant or, 
equivalently, the bag pressure finite and satisfying all other necessary requirements (see 
section 8 above). The key elements of this program were the necessary subtractions at 
all levels. Moreover, one of its attractive features, as emphasized above, is that it is 
the energy density of the purely transversal virtual gluon field configurations which are 
not only stable (no imaginary part), but are being in the stationary state as well, i.e., 
in the state with the minimum of energy (see Fig. 4). That is why it makes sense to 
discuss the "releasing" of the bag constant from the vacuum, more precisely the bag 
energy (B.4). 

From the quantum statistical mechanics point of view, the energy is nothing but 
the pressure multiplied by the volume V in the infinite- volume limit [87] . So the vacuum 
energy E vac in terms of the bag constant is and in GeV units it diverges as follows: 

V o 

E vac = -B YM V = -E YM — o ~ -A 3 , A ->• oo, (D.l) 

cur 

since V j cm 3 ~ A 3 always when the dimensionless UV cutoff A goes to infinity. Evidently, 
in deriving Eq. (D.l) we use the general relation Eym — Bym cm 3 , which is valid in 
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any units for energy (see appendix B above). 

Let us imagine now that we can release the finite portion Eym (B.4) from the 
vacuum in k different places (different "vacuum energy releasing facilities" (VERF)). It 
can be done by n m times in each place, where m = 1, 2, 3...k. Then the releasing energy 
E r becomes 

k 

E r = E YM ^n m . (D.2) 

m=l 

The ideal case (which, however, will never be achieved) is when we could extract a finite 
portion of the energy an infinite number of times and in an infinite number of places. 
So the releasing energy (D.2) might be divergent as follows: 

k 

E r = E Y m x h m / n m ~ A 2 , A — » oo, (D.3) 

m=l 

since the sum over m diverges quadratically in the A — > oo limit, and k ~ A, n m ~ A 
in this case. The difference between the vacuum energy (D.l) and the releasing energy 
(D.3) which is nothing but the remaining in the vacuum energy Er becomes 

E R = E vac -E r = E vac [l + 0(l/\)}, A^oo, (D.4) 

i.e., the QCD vacuum is an infinite and permanent reservoir of energy. The situation is 
even "better" if one takes into account the PT contributions to the vacuum energy (in 
this case the convergence becomes of the order 0(1/A 2 ) in Eq. (D.4), see our preliminary 
work in Ref. [71]). 

That's the vacuum energy is badly divergent is not a mathematical problem. This 
reflects an universal reality. Vacuum is everywhere and it always exists. Quite possible 
that our Universe in general and our real word in particular is only its special type of 
excitation due to the Big Bang. As underlined above, the vacuum is an infinite and hence 
a permanent source of energy. The only problem is how to release the finite portion - 
the bag energy (B.4) and whether it will be profitable or not by introducing some type 
of cyclic process. However, due to huge difference between the two numbers (B.4) and 
(B.5) such a cyclically profitable process may be realistic. "Perpetuum mobile" does not 
exist, but "perpetuum source" of energy does exist, and it is the QCD ground state. 
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